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Abstract 

When a strong electric field is applied to a colloidal suspension, it may cause 
an aggregation of the suspended particles in response to the field. In the case 
of a rotating field, the electrorotation (ER) spectrum can be modified further 
due to the local field effects arising from the many-particle system. To capture 
the local field effect, we invoke the Maxwell-Garnett approximation for the 
dielectric response. The hydrodynamic interactions between the suspended 
particles can also modify the spin friction, which is a key to determine the 
angular velocity of ER. By invoking the spectral representation approach, we 
derive the analytic expressions for the characteristic frequency at which the 
maximum angular velocity of ER occurs. From the numerical caculation, we 
find that there exist two sub-dispersions in the ER spectrum. However, the 
two characteristic frequencies are so close that the two peaks actually overlap 
and become a single broad peak. We report a detailed investigation of the 
dependence of the characteristic frequency and the dispersion strength of ER 
on various material parameters. 

PACS Number(s): 82.70.-y, 87.22.Bt, 77.22.Gm, 77.84.Nh 
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I. INTRODUCTION 



AC electrokinetic pheonomena, such as dielectrophoresis and electrorotation (ER) have 
been investigated to yield the selective manipulation and characterization of biological cells. 
Under the action of external fields, colloidal particles or biological cells in suspensions exhibit 
rich fluid-dynamic behavior as well as dielectric response. It is also interesting to investigate 
their frequency-dependent responses to ac fields, which leads to valuable information on the 
structural (Maxwell- Wagner) polarization effects The polarization is characterized by 
a variety of characteristic frequency- dependent changes known as dielectric dispersion. 

In the last two decades, various experimental tools have been developed to analyze the 
polarization of biological cells - dielectric spectroscopy |§ , dielectrophoresis [§]] and ER f| 
techniques. Among these techniques, conventional dielectrophoresis and ER are usually 
applied to analyze the frequency dependence of translations and rotations of cells in a 
rotating electric field, respectively Moreover, one is able to monitor the cell movements 
with the aid of automated video analysis || as well as light scattering methods 0. In ER, 
a dipole moment is induced in response to the rotating field. Then any dispersion process is 
able to cause a phase shift between the induced dipole moment and the external field vector, 
resulting in a desired torque which causes the cells to rotate. 

In the dilute limit, the ER of individual cell can be predicted by ignoring the mutual 
interaction between cells. However, the cells may be aggregated under the influence of the 
external field. In this case the Brownian motion can be neglected, and thus the system 
becomes non-dilute even though it is initially dilute. Then the role of mutual interaction 
will become dominant. As an initial model, we studied the ER of two approaching spherical 
particles in the presence of a rotating electric field in a recent paper . We showed that when 
the two particles approach and finally touch, the mutual polarization interaction between 
the particles leads to a change in the dipole moment of individual particles and hence the 
ER spectrum, as compared to that of isolated particles. The mutual polarization effects 
were captured via the multiple image method ||. From the results, we found that the 
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mutual polarization effects can lower the characteristic frequency at which the maximum 
ER angular velocity occurs. 

When the volume fraction of the suspension becomes large, the mutual interactions be- 
tween the suspended particles can also modify the spin friction, which is a key to determine 
the angular velocity of ER. For ER of two particles in a rotating electric field, we have suc- 
cessfully applied the spectral representation to calculate the dispersion frequency. However, 
the determination of the spin friction was lacking. 

For two particles, the basic tool is the reflection method [[J, being analogous to the 
multiple image method in electrostatics, but being valid for two particles only. For more 
than two particles, we need a first-principles method, e.g., the Green's function (Oseen 
tensor) formulation. For a dilute suspension, however, one may adopt the effective medium 



theories [|1(J to capture the effective viscosity of a suspension. Thus, in this case, the ER 
spectrum can be modified further. 

Regarding the spectral representation approach |IT| , it is a rigorous mathematical for- 
malism of the effective dielectric constant of a two-phase composite material. It offers the 
advantage of the separation of material parameters (namely the dielectric constant and con- 
ductivity) from the particle structure information, thus simplifying the study. In the present 
work, we will derive the analytic expressions for the characteristic frequencies by using the 
spectral representation approach. 

II. FORMALISM 

We consider homogeneous, spherical particles exposed to a rotating electric field of fre- 
quency /, in which particles of complex dielectric constant Z\ — ei + <Ji/(i27rf) are dispersed 
in a suspension of I2 = ^2 + ^2/ (i^f), with i = Here, e» and <Tj represent the dielectric 

constant and conductivity, respectively. 

For the effective dielectric constant e e of the whole system, we have: 

e e — £2 ?i — £2 , , 

P^— — , (1) 



e e + 2e 2 h + 2e 2 



where p is the volume fraction of partices. Note that the p = limit throughout the paper 
is just the isolated spherical result 0. In this case, the dipole factor of the particle is given 
by 

6 = i^. (2) 

In order to obtain the analytic expression for the characteristic frequency, we resort to 
the spectral representation approach. By introducing a material parameter 
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Eq.(||) admits the form: 
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where 



p = (i-p)(Vx+p-vp) F = (i^pKVg+p + vg 

1 " 6^/WTp ' 2 6y/WTp 
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s 1 = - (2+p-^p(8 + p)j , s 2 = - (2+p+^/p(8 + p)j . 

After some simple manipulations, Eq.(||]) becomes 1?1JT211 

u F i Ae i F 2 Ae 2 

S-Si I + if /td S-S 2 l + lf/fc2 

where the dispersion magnitudes Ae and characteristic frequencies f c being given by 

^ s — t . ^ s — t 

Aei = F x - -, Ae 2 = F 2 - 



(t — Si)(s — SiY ' {t - s 2 ) {s - s 2 ) ' 

1 cr 2 s(* - si) 1 o- 2 s(t - s 2 ) 
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2tt e 2 i(s - si) ' 27r e 2 i(s - s 2 ) ' 

and s = 1/(1 — ei/e 2 ) and £ = 1/(1 — 0i/er 2 ) are the dielectric and conductivity contrasts, 
respectively. Note that we predict two characteristic frequencies f c ± and / c2 . 

Regarding the dynamic viscosity, the noninteracting spin friction expression for an iso- 
lated spherical particle must be modified. For many particles interacting in a suspension, we 
must consider the suspension hydrodynamics. The spin friction is calculated by considering 



the hydrodynamics of a rotating particle in a homogeneous suspension with an effective 
viscosity r\ e . The effective viscosity for a suspension may be obtained from an analogous 
Maxwell-Garnett approximation (MGA) [ |TC|j : 

*^ = P ^, (5) 
Ve + %i) 2 r/i + fr/2 

where rji and 772 are respectively the viscocity of particle and host, and we may take 771 — > 00 
for hard spheres. The dilute limit expression for the effective viscocity is: 

Ve = r}2 + (^p) V2 + (Jp 2 ) rj 2 H . (6) 

This equation is a suitable method to determine the effective viscosity at low concentra- 
tion as it predicts no percolation threshold. The effective medium approximation (EMA) is 
also a useful method to evaluate the effective viscosity of a suspension, which admits the 



form [pL0|| : 



We should remark that a percolation threshold p c = 0.4 is predicted by Eq.(|7|). From this 
equation, the dilute limit expression is 



Ve = V2+ [^Pj V2 + ( 2? ) ^2 H ■ (8) 

To one's interest, Eq.(^) and Eq.(H) attain the same expansion up to first order in p. In 
addition, we should also remark that both Eq.([5D and Eq.(|7|) are valid only at low Reynolds 
number. 

It is known that the rotational angular velocity of a particle is related to the dipole factor 
as follows 

m = - e j^ Im ib] (9) 

where E is the strength of the applied electric field, and Im[- ■ •] means taking the imaginary 
part of [•■■]. 



III. NUMERICAL RESULTS 



During the numerical calculation, set E = 1.0 x 10 A V/m, i] 2 = 1.0 x 10~ 3 Kg/m- s. Note 
e denotes the dielectric constant of free space. 

Fig. 1(a) is plotted versus frequency / for three different e\ at <j 2 = 2.8 x 10~ A S/m, 
e 2 = 80e , p = 0.05 and <j\ = 2.8 x 10~ 2 S/m. Given e\, there is always one peak for angular 
velocity as frequency increases. The frequency located by the peak is just the characteristic 
frequency. We should remark here that two peaks we predict above are located so close that 
only one peak behavior is shown. In our previous work pj] , a pair of particles is investigated 
by considering the mutual depolarization effect, and we found that one peak is dominant 
while the other has a lower value and is located in a quite different position. This shows that 
the particle aggregation effect may enhance the second peak value, and be relocated close 
to the dominant one. Regarding the related quantitative explanation, please refer to Fig. 3 
arising from the spectral representation approach. As e% increases, both the characteristic 
frequency and the angular-velocity peak value decreases. 

Fig. 1(b) is plotted versus frequency / for three different e 2 at a 2 = 2.8 x 10~ 4 S/m, 
p = 0.05, ei = 10eo and o\ = 2.8 x 10 -2 S/m. It is evident that increasing e 2 yields 
decreasing characteristic frequency, but increasing angular-velocity peak value. 

Fig. 1(c) is plotted versus frequency / for three different o\ at a 2 = 2.8 x l0~ A S/m, 
e 2 = 80eo, p = 0.05 and e\ = 10eo- It is evident that both the characteristic frequency and 
the angular- velocity peak value increase for increasing <j\. 

Fig. 1(d) is plotted versus frequency / for three different <j 2 at e 2 = 80eo, p = 0.05 and 
€i = 10eo and o\ = 2.8 x 10~ 2 5/m. For an increasing cr 2 , the characteristic frequency 
increases whereas the angular-velocity peak value decreases concomitantly. 

In Fig. 2, the effect of volume fraction on the effective viscosity is investigated. Clearly, 
at a larger volume fraction, a larger effective viscosity is predicted both for MGA as well as 
EMA. It is evident that, for a given p, EMA predicts a larger effective viscosity than MGA. 

In Fig. 3, we investigate the dependence of the spectral parameters, namely the residues 
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(Fi and F 2 ), the poles (si and s 2 ), as well as the dispersion parameters (Aei and Ae 2 ) on 
the volume fraction p (for MGA only). 

In Fig. 3 (a) and 3(b), as the volume fraction p increases, both F\ and F 2 increase. In- 
creasing volume fraction p leads to decreasing si, but increasing s 2 . Note that the residues 
satisfy the relation Fi + F 2 = —(1 — p)/3. In the two figures, the distinction between Fx and 
i*2, or si and s 2 is quite small at low volume fraction (e.g., p = 0.05, 0.1, and 0.15). 

In Fig.3(c), as the volume fraction increases, Aei first increases, and then decreases, 
while Ae 2 decreases monotonically. Note Ae 2 becomes very close to Aei at low concentration 
region. 

In Fig. 3(d), as the volume fraction increases, f c \ decreases, while / c2 increases concomi- 
tantly. Moreover, it is evident that both characteristic frequencies collapse at low concen- 
tration region. 

From Fig.3, it is readily concluded that at low volume fraction, the angular-velocity 
peaks related to two characteristic frequencies are so close that they can almost be seen as 
overlapped. This is just the reason why only one rotation peak occurs in our figures. We 
believe that such distinction between the two peaks will become recognizable for a high- 
concentration case. 

Fig. 4 is plotted for three different volume fraction p at <j 2 = 2.8 x lO^^S/m, e 2 = 80eo, 
and ei = 10e and a x = 2.8 x 10~ 2 S/m. In this figure, we investigate the effect of particle 
volume fraction on the angular velocity. Increasing volume fraction is able to reduce the 
characteristic frequency as well as the angular-velocity peak value. 

IV. DISCUSSION AND CONCLUSION 

Here a few comments on our results are in order. We considered homogeneous, spherical 
particles exposed to a rotating electric field. As the strength of the rotating electric field 
increases, the particles in a suspension may be aggregated into sheet-like structures. Within 
these structures, the concentration of particles is high. In this case, the local-field effect 
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arising from the many-particle system must be taken into account. In this regard, we 
believe the effective medium theories including both the Maxwell- Garnett approximation 
(also known as the Clausis-Mossotti approximation) and the effective medium approximation 
(also known as the Bruggeman approximation) are good candidate theories. The Maxwell- 
Garnett approximation is well-known to be non-symmetrical and May thus be suitable for 
low concentration. For a higher concentration of particles, we had better use the effective 
medium approximation (EMA). Both methods are valid at low concentrations whereas EMA 
predicts a percolation threshold [llOfl . 

We believe our theory will be valid for the case of low Reynolds numbers. In our present 
study, both the radii and the angular velocities of the particles are small. Thus, we may 
safely use the effective medium theories for the viscosity, which are valid when the particles 
are at rest. 

We realize the appearance of a circular medium flow because all the particles rotate in 
the same direction. The macroscopic spin rate may drive the suspending liquid, leading 
to a decrease of the apparent viscosity of the suspension [13]. Finally, we can apply our 
formalism to single or multi shell objects, like biological cells [12]. 

In summary, we have investigated the electrorotation spectrum of many interacting par- 
ticles including their electric and hydrodynamic interactions. The results showed that the 
electrorotation spectrum can be modified due to the local field effects arising from the many- 
particle system. The hydrodynamic interactions between the suspended particles can also 
modify the spin friction, which is a key to determine the angular velocity of electrorotation. 
By invoking the spectral representation approach, we derived the analytic expressions for the 
characteristic frequency at which the maximum angular velocity of electrorotation occurs. 
We reported a detailed investigation of the dependence of the characteristic frequency and 
the dispersion strength on various material parameters. 
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FIGURES 

FIG. 1. Dependence of rotational angular velocity on ei, e 2 , o\ and 02 versus frequency (for 
the MGA only), (a) cr 2 = 2.8 x W~ 4 S/m, e 2 = 80e , p = 0.05 and <7i = 2.8 x 10 _2 5/m. (b) 
0-2 = 2.8 x 10 _4 5/m, p = 0.05 and ei = 10e and o x = 2.8 x 10~ 2 S/m. (c) a 2 = 2.8 x 10 -4 S/m, 
e 2 = 80e , p = 0.05 and e\ = 10e . (d) e 2 = 80e , p = 0.05 and t\ = 10e and <j\ = 2.8 x 10 _2 S'/m. 

FIG. 2. The MGA and EM A effective viscosity of suspension plotted versus volume fraction p. 

FIG. 3. The spectral parameters Fi, F 2 , s\ and S2, as well as the dispersion parameters, Aei 
and Ae2, plotted versus volume fraction p (for the MGA only). 

FIG. 4. The rotational angular velocity plotted versus frequency for different volume fraction 
p. Parameters used in the plot are e\ = 10eo, ci = 2.8 x 10~ 2 S/m, a 2 = 2.8 x 10 _4 S'/m, e 2 = 80eo 
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